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ABSTRACT 


The  direct  assessment  of  heat  flux  from  the  body  is  a  basic  measurement  in  thermal 
physiology.  Heat  flux  transducers  (HFTs)  are  being  used  increasingly  for  that  purpose  under 
different  environmental  conditions.  However,  questions  have  been  raised  regarding  the  accuracy 
of  the  manufacturer's  constant  of  calibration,  and  also  about  the  effect  of  the  thermal  resistance  of 
the  device  on  the  true  thermal  flux  from  the  skin.  Two  different  types  of  waterproofed  HFTs  were 
checked  for  their  calibration  using  the  Rapid-k  thermal  conductivity  instrument.  A  detailed 
description  of  the  methodology  used  during  the  calibration  is  given.  The  mean  differences  between 
our  calibration  constants  and  the  manufacturers'  constants  were  +20.2  ±  7.1%  (n  «  15)  for 
Thermonetics  Corporation's  HFTs  (San  Diego,  CA)  and  -0.7  +  4.8%  (n  »  12)  for  Concept 
Engineering's  HFTs  (Old  Saybrook,  CT).  The  highly  significant  statistical  difference  in  the  error 
of  calibration  between  the  two  manufacturers  (p  <  0.001)  becomes  an  important  criterion  for  the 
selection  of  HFTfc. 

A  model  capable  of  simulating  a  large  range  of  "tissue"  insulation  was  used  to  study  the 
effect  of  the  underlying  tissue  insulation  on  the  relative  error  in  thermal  flux  due  to  the  thermal 
resistance  of  the  HFIfc.  The  data  show  that  the  deviation  from  the  true  value  of  thermal  flux 
increases  with  the  reciprocal  of  the  underlying  tissue  insulation  (r  *  0.99,  p  <  0.001).  The 
underestimation  of  the  heat  flux  through  the  skin  measured  by  an  HFT  is  minimal  when  the  device 
is  used  on  vasoconstricted  skin  in  cool  subjects  (3  to  13%  error),  but  becomes  important  when 
used  on  warm  vasodilated  subjects  (29  to  35%  error),  and  even  more  important  on  metallic  skin 
mannequins  (  >  60%  error).  In  order  to  optimize  the  accuracy  of  the  heat  flux  measurements  by 
HFTs,  it  is  important  to  recalibrate  the  HFTs  from  Thermonetics  Corporation,  and  to  correct  the 
heat  flux  values  for  the  thermal  resistance  of  the  HFT  when  used  on  vasodilated  tissues. 
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1 .  INTRODUCTION 

The  direct  assessment  of  heat  flux  from  the  body  is  a  basic  measurement  in  thermal 
physiology.  Heat  flux  transducers  (HFTs)  are  being  used  increasingly  for  this  purpose  in 
preference  to  calorimeters  because  of  their  reduced  cost,  their  ease  of  application  and  their 
capability  to  be  used  under  different  and  often  more  realistic  environmental  conditions. 

A  HFT  consists  of  a  sensitive  thermopile  composed  of  many  fine-gauge  thermocouples 
connected  in  series  on  opposite  sides  of  a  flat  matrix  with  a  stable  thermal  resistance  (see  Fig.  1). 
This  thermal  resistance  RT  creates  a  temperature  difference  AT  across  the  matrix  in  the  presence  of 
heat  flux  A 


AT-Rt*H.  (1) 

The  temperature  difference  is  detected  by  the  thermocouple  junctions,  and  the  voltage 
generated  across  the  thermopile  is  proportional  to  the  heat  flux  through  the  HFT.  The  sensitivity  of 
the  device  (mV*m2*W'])  is  proportional  to  the  number  of  thermocouple  junctions  in  series  and  to 
the  thermal  resistance  RT  of  the  matrix. 

Despite  the  simplicity  of  the  design,  questions  have  been  raised  regarding  the  reliability  and 
accuracy  of  HFTs.  Nuckols  and  Piantadosi  (7)  developed  a  calibration  technique  for  HFTs  and 
showed  an  effect  of  temperature  on  the  performance  of  the  HFTs,  as  well  as  significant  differences 
of  up  to  20%  between  the  factory-supplied  calibration  constants  (using  HFTs  from  Thermonetics 
Corporation)  and  those  determined  experimentally.  The  authors  concluded  that  potentially  large 
errors  in  calculated  convective  heat  flow  may  arise  if  transducers  are  not  accurately  calibrated. 
They  suggested  that  HFTs  should  be  carefully  recalibrated  to  avoid  this  error. 

Recently,  Sowood  (12),  using  a  similar  method  of  calibration,  showed  that  the 
experimentally  determined  calibration  constants  of  HFTs  were  the  same  as  those  of  the 
manufacturer  (using  HFTs  from  Concept  Engineering)  minus  approximately  6  Wntf^rniV*1.  He 
concluded  that  the  differences  between  the  two  sets  of  values  were  negligible  and  that  the  devices 
could  be  used  without  recalibration. 

These  two  studies  raise  questions  regarding  a  potential  difference  in  the  accuracy  of  HFT 
calibrations  between  manufacturers.  If  such  a  difference  is  demonstrable  when  using  a  single 
reliable  method  of  calibration  for  both  types  of  HFTs,  the  implication  could  be  important  in  the 
selection  of  HFTs  by  those  who  do  not  have  access  to  calibration  facilities. 

Poppendick  (8),  Gin  et  al.  (6),  and  later  Wissler  and  Ketch  (14)  observed  another  potential 
source  of  error  involved  in  the  use  of  HFTs.  The  design  of  the  HFT  involves  compromise.  The 
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Figure  1:  Temperature  profile  through  a  biological  heat  transfer  system  including  a  heat  flux 
transducer. 
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addition  of  the  thermal  resistance  RT  of  the  HFT  on  the  surface  under  study  will  modify  the 
insulation  of  the  region  and  consequently  change  the  skin  temperature  and  heat  flux.  When  an 
HFT  is  usedron  the  skin  in  water,  where  the  thermal  resistance  between  the  skin  surface  and  the 
environment  is  small,  a  significant  local  increase  in  thermal  resistance  may  occur.  Therefore,  the 
voltage  obtained  from  the  HFT  is  only  representative  of  the  heat  flow  from  the  tissue  immediately 
below  the  disc,  not  that  of  the  adjacent  areas  of  exposed  skin.  Wissler  and  Ketch  (14)  calculated  in 
their  study  that  when  used  in  water,  the  HFT  can  underestimate  the  thermal  flux  by  15  to  20% 
because  of  the  thermal  resistance  of  the  transducer.  In  addition,  Gin  et  al.  (6)  suggested  that  the 
error  could  be  greater  when  HFTs  ore  used  on  vasodilated  skin.  Despite  the  suggestion  that  the 
magnitude  of  the  error  might  be  related  to  the  vascular  state  of  the  underlying  tissue,  no 
information  is  available  in  the  literature  about  the  relationship  between  the  error  in  heat  flux  values 
due  to  the  resistance  of  HFTs  and  the  underlying  tissue  insulation. 

This  study  presents  a  detailed  description  of  the  methodology  used  for  calibrating  HFTs. 
The  recalibration  of  HFTs  from  two  companies  (Thermonetics  Corporation  and  Concept 
Engineering)  is  presented  and  recommendations  are  given  related  to  the  selection  and  use  of  HFTs. 
In  addition,  the  study  gives  new  evidence  for  the  importance  of  correcting  the  thermal  flux  values 
for  the  thermal  resistance  of  the  HFT  when  used  under  certain  conditions.  The  latter  information 
has  been  published  in  condensed  form  in  the  open  literature  (5). 
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2 .  MATERIAL  AND  METHODS 


2.1  Calibration  of  the  Rapid-k  instrument 


2.1.1  The  Rapid*k  instrument 

The  Rapid-k  instrument  (Dynatech  R/D  Company,  Mass)  is  an  apparatus  designed  to 
determine  the  thermal  conductivity  of  materials  in  accordance  with  ASTM  C518  "Specifications  for 
the  Measurement  of  Thermal  Conductivity  by  Means  of  the  Heat  Flow  Meter".  When  calibrated 
properly,  it  can  operate  as  a  reference  instrument  for  the  calibration  of  HFTs. 

Heat  flow  through  a  solid  results  from  having  a  temperature  gradient  in  the  material.  The 
thermal  conductivity  is  a  material  property  which  determines  how  much  heat  flows  through  a  given 
thickness  of  the  material  when  submitted  to  a  temperature  difference.  The  Fourier  linear  heat  flow 
equation  defines  thermal  conductivity  under  steady  state  conditions  as: 


k-JilL  (2) 

AT 

where  k  is  the  thermal  conductivity  of  the  material  in  W^cHnu*1,  li  is  the  heat  flux  through  the 
material  in  W«m'2,  x  is  the  thickness  of  the  material  under  investigation,  and  AT  is  the  temperature 
difference  between  the  two  opposite  surfaces  of  the  material  in  °C 

The  determination  of  thermal  conductivity  using  the  Rapid-k  instrument  is  made  by  placing 
a  sample  of  the  material  (30  cm  wide  x  30  cm  long  x  2  to  5  cm  thick)  between  two  copper  plates 
which  are  maintained  at  known  temperatures.  Because  of  the  temperature  difference  between  the 
surfaces,  heat  will  flow  through  the  sample  from  the  hot  side  to  the  cold  side  (see  Fig.  2A).  The 
quantity  of  heat  flowing  through  the  sample  is  measured  by  a  very  sensitive  heat  flow  meter  located 
between  the  sample  and  the  cold  plate  (the  term  "meter"  is  used  in  reference  to  the  internal  HFT  of 
the  Rapid-k  to  avoid  confusion  with  the  HFT  under  test).  The  heat  flow  through  the  central  10  cm 
x  10  cm  area  of  the  sample,  called  the  effective  area,  is  unidirectional  and  uniform  between  the  hot 
and  the  cold  plate. 

If  the  k  value  of  the  sample  is  known,  the  value  of  Hcan  be  calculated  knowing  x  and  AT. 
The  relationship  between  the  voltage  output  of  the  heat  flow  meter  and  H  can  be  determined  by 
applying  different  ATs  to  the  copper  plates.  This  is  the  basic  principle  of  the  Rapid-k  calibration. 
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.Hotplate 

30  cm _ ,  Standard  specimen 


H 


Figure  2:  Cross-sectional  (A)  and  top  (B)  views  of  the  Rapid-k  test  set-up  and  placement  of  the 

thermocouples  for  the  measurement  of  AT  between  the  two  opposite  plate  surfaces  of  the 
thermal  resistance  standard.  The  thermocouple  junctions  are  fixed  inside  the  effective 
area  where  the  value  A  is  measured  by  the  heat  flow  meter  of  the  Rapid-k.  The  gaps 
between  the  plates  and  the  standard  on  A  are  present  only  for  the  clarity  of  the  figure. 
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To  ensure  optimal  performance  of  the  Rapid-k,  the  heat  flow  meter  of  the  instrument  must 
be  calibrated  against  a  thermal  resistance  standard.  For  that  purpose,  two  thermal  resistance 
standard  specimens  were  borrowed  from  the  Thermal  Performance  Section  of  the  National 
Research  Council  of  Canada  (NRC). 


2.1.2  The  National  Research  Council  thermal  resistance  standards 

Table  1  presents  the  calibration  data  provided  for  the  standards  (MDGB  Transfer  Standards 
#357-172-A,B)  by  NRC.  The  two  calibrated  standard  specimens  had  identical  calibrations: 

R  » (0.88258 1  ±  1 .40447E-03)  -  (2.87934E-03  ±  5.34837E-05)*T  (3) 

where  R  is  the  thermal  resistance  of  the  standard  in  °Om2»W’1,  and  T  is  the  mean  temperature  of 
the  standard  in  °C.  The  thicknesses  of  the  acrylic  framed  standards  were  26.03  mm  for  specimen 
A  and  26.01  mm  for  specimen  B. 


2.1.3  Rapid-k  calibration 

This  calibration  intended  to  establish  the  relationship  between  the  heat  flow  meter  output  of 
the  Rapid-k  (V\  in  mV)  and  the  calculated  heat  flux  (H  in  W*m"2)  using  the  standard  resistance 
specimen  of  known  thermal  resistance  (R  in  °C*m2»W'1).  The  following  equations  were  used  io 
calculate  the  heat  flux  through  the  standard  when  exposed  to  a  temperature  difference  AT: 

ft  =  R  *  AT  (4) 


or 


ft  =  (0. 88258 1  -2.87934E-03»T)*(Thot-Tcold)  (5) 

where  T  is  the  mean  temperature  of  the  standard.  The  temperature  difference  AT  was  measured 
accurately  with  six  fine  gauge  thermocouples  (AWG  40)  calibrated  against  a  quartz  thermometer 
(HP  2804A)  to  an  accuracy  of  0.07  °C  in  a  well-stirred  temperature  controlled  water  bath  between  5 
and  80°C.  Preliminary  experiments  have  shown  that  by  using  the  thermocouples  embedded  in  the 
copper  plates  of  the  Rapid-k  to  measure  the  AT,  an  error  of  approximately  4%  is  introduced  in  the 
calculation  of  ft.  Three  thermocouples  were  fixed  on  each  surface  of  the  standard  within  the  10  cm 
x  10  cm  effective  area  where  the  heat  flow  meter  of  the  Rapid-k  is  located  (see  Fig.  2B).  The  mean 
temperature  of  the  specimen  was  calculated  from  the  mean  value  of  the  six  thermocouple  readings. 


Table  1:  Thermal  resistance  data  provided  by  the  National  Research  Council  for  the  MDGF 
Transfer  Standards  #357- 1 72- A, B . 


Mean  temperature 
of  standard 

(*C) 

Measured  thermal 
resistance 

fC-mMV-1) 

Calculated  thermal 
resistance 

Difference  between  measured 
and  calculated  values 

24.17 

.8130 

.812987 

1.31726E-05 

20.14 

.8283 

.824501 

3.60942E-03 

23.72 

.8185 

.814282 

4.21751E-03 

23.93 

.8129 

.813678 

-7.77900E-04 

0.90 

.8788 

.879989 

-1.18911E-03 

13.07 

.8450 

.844948 

5.24521E-05 

20.53 

.8223 

-  .823468 

-1.16766E-03 

22.27 

.8167 

.818458 

-1.75756E-03 

24.26 

.8128 

.812728 

7.23600E-05 

26.10 

.8071 

.807171 

-7.05710E-05 

28.01 

.8001 

.801930 

-1.83016E-03 

37.31 

.7736 

.775152 

-1.55222E-03 

47.95 

.7449 

.744518 

3.83073E-04 
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The  calibration  of  the  Rapid-k  was  performed  for  values  of  H  ranging  from  0  to  100  W»m'2. 
The  upper  extremity  of  the  range  was  limited  by  the  low  melting  point  temperature  of  the  acrylic 
"  used  to  frame  the  standard  specimen.  Before  each  calibration  test,  the  specimen  was  dried  in  an 
oven  for  15  hours  at  70°C,  since  the  R  value  given  by  NRC  was  for  a  dry  specimen.  The  values  of 
temperature  and  voltage  were  considered  to  be  at  thermal  stability  when  the  ratio  of  voltage  from 
the  heat  flow  meter  over  AT  varied  less  than  1  %  over  a  period  of  20  minutes. 

Figure  3  depicts  the  calibration  curve  of  the  Rapid-k,  representing  the  relationship  between 
the  calculated  A  in  W*m-2  and  the  voltage  output  of  the  heat  flow  transducer  in  mV  The  highly 
significant  (p  <  0.001)  linear  relationship  can  be  described  as  follows: 

ft  «  27.5092*  V-0.0044  (6) 

r  *  0.999 

By  varying  the  temperature  of  the  heat  flow  meter  of  the  Rapid-k  from  5  to  30°C  (i.e.,  the 
cold  plate  temperature),  a  variation  in  the  voltage  output  of  less  than  1%  was  observed  for  the  same 
calculated  A  values.  It  can  therefore  be  concluded  that  the  relationship  between  the  calculated  A 
and  V  is  temperature  independent  at  least  for  this  calibration  temperature  range. 


2.2  Characterization  of  HFTs 

Two  waterproof  versions  of  HFfls  available  on  the  market  (one  made  by  Thermonetics 
Corporation  and  the  other  by  Concept  Engineering)  were  chosen  for  evaluation  because  they  offer 
the  best  potential  for  the  measurement  of  heat  flow  in  a  variety  of  environments. 

Figures  4  and  5  depict  the  characteristics  of  the  two  types  of  HFT  discs.  Tne  HFTs  from 
Thermonetics  are  1.1  to  2.1  times  more  sensitive,  but  more  fragile  than  those  from  Concept.  We 
observed  that  the  Thermonetics  HFT  is  susceptible  to  mechanical  breakdown  at  the  junction 
between  the  disc  and  the  wire.  Because  the  Concept  HFT  is  provided  without  the  wire  (the  wire  is 
optional),  the  strength  of  the  junction  depends  upon  the  material  used  to  seal  the  connection 
between  the  HFT  and  the  wire.  Another  difference  between  the  two  HFT  discs  is  the  location  of 
the  thermistor  on  the  disc.  The  placement  of  the  thermistor  on  the  edge  of  the  Thermonetics  HFT 
makes  it  more  susceptible  to  thermal  measurement  error  due  to  the  edge  effect  described  by  Wissler 
and  Ketch  (14).  In  both  cases  however,  the  thermistor  measures  the  temperature  of  the  skin  under 
the  HFT,  which  can  be  different  from  the  temperature  of  the  skin  not  covered  by  the  HFT  (14). 


Voltage  output  (mV) 


Figure  3: 


Calibration  curve  of  the  Rapid-k  instrument  representing  the  voltage  output  from 
flux  meter  of  the  Rapid-k  (mV),  and  the  calculated  heat  transfer  H  (W*m',fc) . 
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Model  #:  HA13-18-10*P(C) 

Factory  calibration  constant:  60  to  100 
Sensitivity:  10  to  17 
Accuracy:  10% 

Impedance:  162  ohms 

Factory  stated  thermal  conductivity:  0.225  W’m”1,  *  C-1 

Thermistor  type:  YSI  44018  (temperature  range:  -30  to  105* C) 

price:  US  $170. 00/unit 

Time  delay  to  order:  several  months 

Address  of  manufacturer:  Thermonetics  Corporation 

Box  9112  San  Diego, 

California  92109 
U.S.A. 

tel:  (619)  488-2242 


^  Thermistor 


Figure  4:  Characteristics  of  the  heat  flux  transducer  from  Thermonetics  Corporation. 


Thermistor 


Model  #:  FR-040-TH44018  (without  wire) 

Factory  calibration  constant:  11C  to  125  W,m-2,mV~1 
Sensitivity:  8  to  9  f.tV'in2-W~x 
Time  constant:  7  sec 
Impedance:  175  ohms,  nominal 

Factory  stated  thermal  conductivity:  0.4615  W'in-1’  *  C-1 

Thermistor  type:  YSI  44018  (temperature  range:  -30  to  105  *  C) 

Price:  US  $190.00/unit 

Time  delay  to  order:  several  weeks 

Address  of  manufacturer:  Concept  Engineering 

43  Ragged  Rock  Road 
Old  Saybrook 
CT  06475  U.S.A. 
tel:  (203)  388-5560 


Figure  5:  Characteristics  of  the  heat  flux  transducer  from  Concept  Engineering. 


2.3 


Calibration  of  HFTt 


r  The  method  of  calibration  for  the  HFTs  presented  in  this  section  is  an  adaptation  of  the 
method  described  by  Nuckols  and  Piantadosi  (7). 

During  the  calibration  tests,  the  HFR  were  held  firmly  within  a  test  bed  in  the  test  section 
of  the  Rapid-k,  The  reasons  for  using  a  test  bed  were:  1)  to  create  a  range  of  heat  flow  through  the 
HFTs  similar  to  that  expected  from  the  body  under  different  environmental  conditions  (from  0  to 
500  W«m‘2);  2)  to  optimize  the  homogeneity  of  heat  flow  around  the  HFTs;  and  3)  to  ensure  good 
surface  contact  with  the  HFTs.  The  criteria  for  the  selection  of  the  test  bed  materials  were:  1) 
resistance  to  heat  deformadon;  2)  a  value  of  thermal  conductivity  close  to  the  k  value  of  the  HFTs 
(between  0.2  and  0.4  W^C’nn'1);  and  3)  some  practical  considerations  (availability,  cost  and  ease 
of  machining).  Two  different  materials  met  those  requirements:  Teflon,  with  a  k  value  of  0.24 
and  hardboard  wood  with  k  «  0.20  (1).  Figure  6  shows  the  test  bed  used 

for  the  calibration  tests.  It  was  comprised  of  two  parts;  namely,  the  test  bed  matrix  and  the  HFT 
support.  Because  of  the  cost  of  teflon,  only  the  HFT  support  was  made  of  this  material.  Four 
HFTs  were  positioned  in  the  HFT  support  within  the  central  10  cm  x  10  cm  effective  area  where 
is  measured  by  the  heat  flow  meter  of  the  Rapid-k.  They  were  held  in  place  by  a  layer  of  the 
same  type  of  tape  as  that  used  to  fix  the  HFTs  over  the  skin  (Blenderm  surgical  tape,  3M 
Corporation).  Tb  optimize  the  heat  transfer  from  the  test  bed  to  the  HFTs,  the  interfaces  were  filled 
with  a  thermally  conductive  grease  (Thermal  compound  part  no.  120-2,  Wakefield  Engineering, 
Wakefield,  MA).  Figure  7  presents  the  complete  test  setup  of  rhe  system  within  the  effective  area 
of  the  test  section  of  the  Rapid-k  as  used  during  a  calibration  run. 

During  the  calibrations,  the  transducers  were  exposed  to  heat  flows  ranging  from  0  to  500 
W*m'2  by  varying  the  temperature  of  the  hot  plate  from  30  to  95°C.  The  temperature  of  the  cold 
plate  was  maintained  constant  at  approximately  30°C  to  minimize  the  effect  of  temperature  on  HFT 
responses  and  to  avoid  the  change  in  k  value  of  the  HFTs  with  temperature.  The  HFT, 
thermocouples  and  Rapid-k  heat  flow  meter  output  voltages  were  read  with  a  computer  controlled 
data  acquisition  system  (HP-85  computer,  HP-3495  scanner,  HP-3455A  digital  voltmeter; 
Hewlett-Packard  Corporation).  After  each  minute  of  scanning,  a  mean  value  was  calculated  for 
each  transducer.  The  program  used  for  the  calibration  of  the  HFTs  is  presented  in  Annex  A.  Only 
values  obtained  during  steady  state  were  used  to  calculate  the  HFT  calibration  constants.  Thermal 
steady  state  was  considered  established  when  the  ratio  of  measured  heat  flux  from  the  Rapid-k  (H) 
to  temperature  differential  (AT  -  Thot  -  Tcold)  did  not  change  by  more  than  1%  over  a  20-min 
period.  The  calculations  of  the  HFT  calibration  constants  were  made  by  transforming  the  voltage 
output  of  the  Rapid-k  heat  flow  meter  into  heat  flow  H  using  the  previously  established  calibration 
curve  of  the  Rapid-k,  and  by  using  the  voltage  output  of  the  HFT  (mV)  under  investigation  as 
follows: 


Test  bed  matrix 
(hardboard) 


Effective  area 

HFT  support 
(teflon) 


Figure  6:  Representation  of  the  test  bed  used  for  the  calibration  of  the  HFTs.  It  was  comprised  of 
the  test  bed  matrix  and  the  HFT  support.  Four  HFTs  were  positioned  on  the  HFT 
support  in  the  effective  area. 


Effective  area 

yv. 


Hot  plate 
Test  bed  matrix 

HFT 

Thermal  conductive  grease 

HFT  support 

Tape 

Heat  flow  meter  of  Rapld-k 
Cold  plate 


H 


Figure  7:  Representation  of  the  effective  area  of  the  test  section  for  the  calibration  and  the 

determination  of  k  value  of  the  HFTs.  T ^ot,  T[  and  Tu  represent  the  location  of  the 
thermocouples  used  to  determined  the  HFT  thermal  conductivities. 
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calibration  constant  — ££P-'k 


(7) 


where  the  calibration  constant  is  expressed  in  Wmr^nnV*1. 

lb  define  the  effects  of  temperature  on  HFT  response,  four  transducers  were  subjected  to 
constant  heat  flows  but  to  six  different  temperatures  ranging  from  10°C  to  60°C  by  varying  the  cold 
plate  temperature  along  with  the  hot  plate  temperature.  In  addition,  the  reproducibility  of  the 
calibration  method  was  verified  by  the  recalibration  of  four  HFR  over  a  two-week  period. 

Verification  of  the  consistency  of  our  calibration  method  was  obtained  using  a  double  layer 
HFT  test.  This  test  consisted  of  measuring  the  value  of  heat  flow  passing  through  two 
superimposed  HFTs  sandwiched  between  two  pieces  of  hard  board  (see  Fig.  8).  Since  the  heat 
flow  is  the  same  through  the  two  HFTs,  the  calculated  values  of  ti,  using  the  calibration  constants 
determined  by  the  Rapid-k  and  the  signal  output  of  each  HFT,  should  be  equal  if  the  calibration 
method  is  consistent 


2.4  Determination  of  HFT  thermal  conductivity 

The  determination  of  the  thermal  conductivity  value  for  the  HFTs  was  performed  at  the 
same  time  as  the  calibration  tests  (see  Fig.  7).  For  each  HFT  in  the  test  bed,  two  calibrated  fine 
gauge  thermocouples  (AWG  40)  were  fixed  on  the  opposite  surfaces  of  the  transducer  (Tt  and  Tu). 
The  k  values  of  the  HFTs  were  calculated  using  the  Fourier  linear  heat  flow  equation  as  follows: 

k»-...&ga —  (8) 

Tr  Tu 

where  k  is  the  thermal  conductivity  of  the  HFT  and  the  tape  in  W*°C*  m*1  is  the  thickness  of  the 
HFT  plus  the  tape  in  m,  and  Tt  -  Tu  is  the  temperature  difference  between  the  two  opposite 
surfaces  of  the  HFT  including  the  tape  in  °C.  The  tape  is  included  in  the  calculation  of  the  k  value 
because,  when  used  over  the  HFT  to  fix  the  transducer  on  the  skin,  it  acts,  like  the  HFT,  as  an 
insulative  material.  The  thermal  resistance  of  the  HFT,  RT,  is  x/k  and  is  in  units  of  °Om2*W'1. 

2.5  Correction  for  the  thermal  resistance  of  HFTs 

2.5.1  Validation  of  the  proposed  equations 
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Hot  plate 
HFT  1 

Thermal  conducting  grease 


HFT  2  - 


Heat  flow  meter  of  Rapid-k 
Cold  plate 


H 


Figure  8:  Representation  of  the  effective  area  of  the  test  section  for  the  double  layers  HFT  tests. 
The  HFT  support  is  made  with  two  1.3  cm  thick  layers  of  hardboard  wood. 
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Three  different  equations  have  been  proposed  in  the  literature  to  correct  the  measured 
values  of  thermal  flux  for  the  effect  of  the  thermal  insulation  of  the  HF'ft  (6,  13,  14).  The 
proposed  equations  are  the  following: 


i)  Strong  etal.  (13) 


H. 


meas 


(T,-Tw)-Rt«H 


meas 


ii)  rearranged  from  the  equation  of  Gin  et  al.  (6) 


(9) 


H. 


(TC-T*)*H 


meas 


VTsk 


corr 


(Tc*-T, ) 


(Tc  ■  Tu)/Hmeas  -  Rt 


iii)  Wissler  and  Ketch  ( 1 4) 


GO) 


H  * _ 322! _ _  (11) 

^  1'Hmeas  -RT/(Tc*Tw) 

where  is  the  thermal  flux  through  the  skin  corrected  for  the  insulating  effect  of  the  HFT 
(W»m*2),  ltmas  is  the  measured  heat  flow  value  (W*m'2),  Tc  is  the  core  temperature  of  the 
segment  for  which  the  heat  flux  is  measured  (°C),  Tsk  is  the  temperature  of  the  uncovered  skin  of 
the  segment  (°C),  Tt  is  the  temperature  of  the  skin  under  the  HFT  (°C),  Tu  is  the  temperature  of  the 
upper  surface  of  the  HFT  (°C),  Tw  is  the  water  temperature  (°C),  and  Rr  is  the  thermal  resistance 
of  the  HFT  (°Om2*W*1) . 

Although  all  the  above  equations  are  theoretically  equivalent  (see  Annex  B  for 
demonstration),  some  of  them  differ  in  practice  due  to  the  difficulties  of  obtaining  accurate  values 
for  the  parameters.  Application  of  the  equations  to  experimental  data  has  shown  that  the  Strong 
equation  (13)  is  too  sensitive  to  temperature  measurement  errors  (as  low  as  0.05°C),  which  makes 
it  impractical.  On  the  other  hand,  the  equations  rearranged  from  Gin  et  al.’s  study  (6)  require 
measurements  of  the  temperatures  of  the  surfaces  of  the  HFTs.  These  are  difficult  to  measure 
accurately,  and  therefore  can  easily  lead  to  error  when  is  calculated.  The  application  of 
Wissler  and  Ketch's  (14)  equation  to  experimental  data  has  shown  it  to  be  useable.  The  only 
practical  limitation  of  this  equation  is  the  measurement  of  the  parameter  Tc,  the  core  temperature  of 
the  segment  under  consideration.  If  it  cannot  be  measured,  one  might  obtain  an  estimate  from 
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mathematical  modelling. 


2.5.2  Effect  of  the  insulation  of  the  underlying  body  tissue  on  the  relative  error 
in  thermal  flux 

The  equation  proposed  for  the  correction  of  the  thermal  resistance  of  HFTs  (i.e.,  equation 
iii)  shows  that  the  correction  factor  in  thermal  flux  is  proportional  to  the  thermal 

resistance  of  the  transducer  (RT),  but  inversely  proportional  to  the  total  thermal  resistance  from 
core  to  ambient  (Tc-Ta  as  expected.  In  addition,  Gin  et  al.  (6),  and  later  Wissler  and 

Ketch  (14),  suggested  that  the  relative  error  in  thermal  flux  is  also  dependent  upon  the  underlying 
insulation  of  the  body  tissue.  When  the  HFTs  are  used  undo*  cold  stress  causing  cutaneous 
vasoconstriction,  the  error  in  thermal  flux  seems  unimportant  (6),  but  when  used  under  conditions 
that  do  not  produce  vasoconstriction,  the  error  in  thermal  flux  may  approach  15%  (14). 

However,  to  our  knowledge,  no  additional  information  is  available  in  the  literature  on  the 
effect  of  the  underlying  insulation  of  the  body  tissue  on  the  relative  error  in  thermal  flux  due  to  the 
thermal  resistance  of  the  HFT. 


2.5.2.1  The  variable-R  model 


In  order  to  assess  the  relationship  between  the  underlying  tissue  insulation  and  the  error  in 
thermal  flux  due  to  the  thermal  resistance  of  HFT,  a  model  system  was  developed  to  create  a  large 
range  of  "tissue"  thermal  resistance  values  (see  Fig.  9).  The  variable-R  model  consisted  of  an 
insulated  water-filled  copper  box  (maintained  at  37°C  with  a  proportional  temperature  controller) 
immersed  in  a  water  bath  maintained  at  30°C  with  a  second  proportional  temperature  controller. 
The  thermal  insulation  of  the  model  was  varied  by  changing  the  material  type  (polystyrene, 
neoprene,  hard  board  wood,  aluminum  and  copper)  and  thickness  0  mm  to  13  mm)  on  one  wall  to 
give  insulation  values  ranging  from  2.7X10"4  to  4.4X10'1  °Om2*W*1.  For  each  setting,  the  thermal 
flux  through  the  system  was  measured  with  3  calibrated  HFTs  (either  from  Thermonetics  or 
Concept  Engineering)  fixed  on  the  external  surface  of  the  insulation  with  Blenderm  surgical  tape. 
The  interfaces  between  the  copper  wall  of  the  variable-R  model  and  the  insulative  material,  and 
between  the  insulative  material  and  the  HFT,  were  filled  with  thermally  conductive  grease.  The 
"skin"  temperature  of  the  model  (Tsk)  was  measured  with  a  fine  type  T  calibrated  thermocouple. 
The  thermal  insulation  of  the  model  ( Rm ,  in  °C»m2*W'1)  was  calculated  as  follows: 


Rm  = 


Tc  *  Tsk 


H 


corr 


(12) 
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(Tv  =  30°C) 


H 


Tape 

HFT 

Film  of  thermally 
conductive  grease 

Insulation  material 
Copper  wall 


Figure  9:  The  variable-R  model  used  for  the  study  of  the  relationship  between  the  underlying 
tissue  insulation  and  the  relative  error  in  thermal  flux  due  to  the  thermal  resistance  of 
HFTfc.  The  thermal  resistance  of  the  model  wall  (copper  box)  was  varied  by  using 
different  insular' ve  material  and  thickness. 
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where  Tc  was  the  core  temperature  of  the  model  (37°C),  was  calculated  by  using  Wissler  and 
Ketch's  (14)  equation.  The  thermal  insulation  values  were  calculated  using  the  mean  values  of 
thermal  flux  and  temperature  recorded  during  the  last  10  minutes  of  each  experiment,  when  thermal 
stability  was  achieved. 


2.S.2.2  Tests  with  human  forearms  immersed  in  water 

The  objective  of  these  tests  was  to  validate  the  data  from  the  variable-R  model  with 
experimental  data  from  immersions  of  human  forearms  exposed  at  different  water  temperatures. 

The  methods  used  during  the  study  were  approved  by  DCIEM  Human  Ethics  Committee. 
The  experiments  involved  ten  healthy  male  subjects  between  18  and  30  years  of  age.  Each 
experiment  consisted  of  immersing  the  forearm  in  a  water  bath  for  3  hours  at  a  constant 
temperature  ranging  between  15  and  36°C.  These  water  temperatures  were  chosen  in  order  to 
create  a  large  range  of  tissue  insulation.  Each  subject  experienced  two  different  temperatures 
chosen  randomly,  one  for  each  forearm,  over  a  two-week  period.  During  the  experiments,  the 
subjects  were  in  a  sitting  position,  and  the  ambient  temperature  of  the  room  was  25  ±  1°C.  A  total 
of  20  experiments  were  performed. 

The  muscle  temperature  (Tm)  was  continuously  monitored  during  the  immersions  at  each  5 
mm  of  tissue  from  the  longitudinal  axis  of  the  forearm  (determined  from  computed  tomography 
scanning)  to  the  skin  on  the  proximal  third  of  the  limb  (superficial  muscle:  flexor  carpi  ulnaris ; 
deep  muscle:  flexor  digitorum  prolundus).  The  temperature  of  the  muscle  was  recorded  with  a 
fine  calibrated  multicouple  probe  implanted  under  local  anaesthesia  (skin  refrigerant,  Ethyl 
Chloride;  Graham-Field,  N.  Y.)  by  using  a  method  detailed  elsewhere  (4).  The  temperature  values 
were  corrected  for  the  thermal  conductivity  effect  along  the  wires  of  the  probe  (4).  The  maximal 
value  of  Tm,  always  found  at  the  longitudinal  axis  of  the  forearm,  was  used  to  represent  Tc  the 
core  temperature  of  the  segment.  The  other  physiological  parameters  recorded  during  the 
immersions  were  the  heat  flux  through  the  skin  of  the  forearm  with  a  recalibrated  Thermonetics 
HFT,  and  the  skin  temperature  with  a  fine  calibrated  40  gauge  thermocouple  probe.  Both  the  HFT 
and  the  thermocouple  were  fixed  a  few  millimetres  beside  the  site  of  the  muscle  temperature 
measurements  using  Blenderm  surgical  tape. 

The  values  were  calculated  by  using  the  Wissler  and  Ketch  correction  previously 
described  (14).  All  thermal  flux  and  temperature  values  used  in  the  calculations  were  from  the  last 
15  minutes  of  each  experiment,  when  thermal  stability  was  achieved. 
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2.6  Statistical  analyses 

The  different  linear  and  curvilinear  regression  analyses  and  independent  t-tests  presented 
throughout  this  study  were  performed  wi*h  the  BMDP  statistical  programs  (BMDP  Statistical 
Software,  1983,  Los  Angeles,  California).  The  values  are  expressed  as  mean  ±  S.D. 
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3.  RESULTS 


3.1  Calibration  of  HFTs 

All  HFTs  from  Thermonetics  end  Concept  showed  linear  responses  over  the  calibration 
heat  flux  range  of  0  to  500  W*m'2.  The  correlation  between  the  HFT  voltage  output  and  the 
Rapid-k  heat  flow  is  highly  significant  (p  <  0.001)  for  each  HFT  calibrated.  Each  transducer 
response  had  a  different  slope,  but  the  linearity  of  the  responses  and  the  nearly  zero  intercepts 
support  die  use  of  a  constant  calibration  factor  for  each  transducer  over  the  test  range  of  heat  flow. 
An  example  of  the  linearity  of  the  response  is  presented  in  Fig.  10  for  4  HFTs  from  Thermonetics. 

The  calibration  constants  determined  by  the  Rapid-k  for  15  HFTs  from  Thermonetics  Corp. 
and  12  HFTs  from  Concept  Engineering,  and  those  constants  provided  by  the  manufacturers  are 
presented  in  Table  2.  The  mean  difference  between  our  calibration  constants  and  the  factory 
constants  is  20.2  ±7.1%  (8.3  to  32.2%)  for  the  Thermonetics  HFTs  and  -0.7  ±  4.8%  (-8.6  to 
8.8%)  for  the  Concept  HFTs  (4.0  ±  2.5%  in  absolute  value  ignoring  the  signs).  The  mean 
difference  value  is  significantly  different  from  0  for  the  Thermonetics  HFTs  (p  <  0.001),  but  not 
for  the  Concept  HFTs  (p  >  0.05).  Both  mean  difference  values  are  significantly  different  from  one 
another  (p  <  0.001).  The  repeat  calibrations  of  4  HFIk  over  a  two-week  period  gave  a  mean 
variability  of  less  than  2%  from  the  first  determination,  which  confirms  the  reproducibility  of  our 
technique.  In  addition,  the  calculated  heat  flow  values  from  two  superimposed  HFTs,  using  the 
respective  calibration  constants  and  voltage  outputs, gave  for  three  tests  a  mean  difference  between 
the  two  values  of  less  than  2%.  These  results  confirmed  the  consistency  of  the  calibration  between 
HFTs.  In  addition,  as  shown  in  Fig.  1 1,  the  temperature  of  the  HFT  had  no  significant  effect  (p  > 
0.05)  on  the  calibration  constant  for  the  4  HFIk  tested  over  a  range  of  temperatures  between  10 
and60°C 


3.2  Thermal  conductivity  of  HFTs 

The  mean  value  of  the  thermal  conductivity  of  the  HFTs  at  30°C  determined  by  the  Rapid-k 
is  0.194  ±  0.018  W»°CInn'1  (n  *  17)  for  the  Thermonetics  HFTs  and  0.316  ±  0.015  W*°C*l*m'1  (n 
»  15)  for  the  Concept  Engineering  HFft.  The  two  values  are  significantly  different  (p  <  0.001) 
and  specific  to  the  construction  characteristics  of  each  type  of  HFT.  The  mean  difference  between 
the  k  value  determined  by  the  Rapid-k  and  the  factory  value  is  16.0%  for  the  Thermonetics  HFTs, 
and  46.0%  for  the  Concept  HFTs . 


H  (W*m 
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Figure  10:  Examples  of  the  linear  relationship  between  the  HFT  voltage  output  and  the  Rapid-k 
heat  flux  for  4  HFTs  from  Thermonetics  Corporation. 


Tfcble  2:  List  of  the  thermal  constants  for  15  KFTs  from  Thermonetics  Corporation  and  12  HFTs 
from  Concept  Engineering,  determined  by  the  manufacturers  and  during  the  recalibration 
tests. 


A.  HFTs  from  Thermonetics  Corporation 


Serial 

number 


Thermal  constant 
from  manufacturer 
(\V*m~2*mV*’1} 

Thermal  constant 
from  recalibration 
(W*m"2imV_1) 

Percentage 
of  deviation 
(%) 

64.7 

55.316 

17.0 

69.4 

56.357 

23.1 

71.0 

62.480 

13.7 

69.4 

53.556 

18.5 

71.0 

57.912 

22.6 

69.4 

57.358 

21.0 

69.i 

82.427 

11.2 

72.5 

66.925 

8.3 

72.5 

62.902 

15.3 

80.4. 

61.535 

30.7 

94.6 

77.515 

22.0 

71.0 

60.962 

16.5 

74.1 

56.458 

31.3 

93.0 

70.332 

32.2 

72.5 

60.900 

19.1 

Mean  20.2:k7.1 

B.  HFTs  from  Concept  Engineering 


Serial 

Thermal  constant 

Thermal  constant 

Percentage 

number 

from  manufacturer 

from  recalibration 

of  deviation 

(Wm“2,mV_l) 

1  (W*m-2*mV-1) 

(%) 

1118 

111.04 

116.758 

-4.9 

1122 

123.66 

120.205 

2.9 

1123 

123.66 

126.817 

•2.5 

1120 

123.66 

121.355 

1.5 

1119 

118.29 

123.095 

-3.9 

JUS 

117.03 

128.000 
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Figure  1 1 :  Effect  of  temperature  on  the  calibration  constant  determined  during  the  recalibration 
tests  for  4  HFTs  from  Thermonetics  Corporation. 
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3.3  Relationship  between  the  underlying  "tissue"  insulation  and  the  relative 
error  in  thermal  flux  due  to  the  thermal  resistance  of  the  HFT 


3.3.1  The  variable«R  model 

Eight  experiments  using  3  HFTs  from  Thermonetics  Corp.  and  7  experiments  using  3 
HFTs  from  Concept  were  performed  using  the  variable-R  model  in  order  to  establish  the 
relationship  between  the  underlying  "tissue"  insulation  and  the  relative  error  in  thermal  flux  due  to 
the  thermal  resistance  of  the  HPT.  Fig.  12  presents  the  relationship  performed  on  mean  values  (3 
HFTs  for  each  data  point)  between  the  ratio  ficcJ^meas'  which  represents  the  correction  factor, 
and  the  ratio  RfiRm  for  all  the  range  of  Rm  studied.  No  significant  difference  (p  >  0.05)  exists 
between  the  relationship  determined  for  each  type  of  HFT.  Tire  curvilinear  relationship  is  best 
described  for  both  types  of  HFTs  by  the  following  equation: 


H  R  R  2  R  3 

*  1.721+1.145  log^l.  +0.688  (log-^ )  +0.156  (log-r^-  )  (13) 

HmCas  Kfn  Km  K"> 

r2  «  0.998 

When  only  the  physiological  range  for  the  ratio  Rj/Rm  is  presented  (RiVRm  <  0.6;  see  ref. 
3),  the  relationship  between  the  ratio  and  the  ratio  Rj^Rm  becomes  linear  (see  Fig. 

13).  The  linear  relationship  is  best  described  for  both  types  of  HFTs  by  the  following  equation: 


corr 

• 

H 

meas 


*  1.004+0.908 


r2  «  0.999 


(14) 


3.3.2  Immersion  of  the  forearm  in  water 

Thble  3  presents-  the  temperature  and  heat  flux  data  for  the  last  15  minutes  of  3  hour 
immersions  of  the  forearm  in  water  at  temperatures  ranging  between  15  and  36°C  (n  «  20).  As 
expected,  the  core  (Tc)  and  the  skin  (Tsk)  temperatures  of  the  forearm  increase  with  the  water 
temperature  (Tw),  and  the  measured  values  of  thermal  flux  from  the  skin  are  inversely 
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proportional  to  Tw.  The  mean  values  of  thermal  flux  corrected  for  the  thermal  resistance  of  the 
HFT  (ticon)  are  9.8  ±  1.6%  higher  than  for  Tw  £  30°C,  and  25.3  ±  7.0%  higher  than 

fimeas  ^  ^ 


Ratio  H 
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5 


4  H 


•  Thermonetics 
o  Concept  Eng. 


Figure  12:  Relationship  between  the  ratio  iloonliimtas,  which  represents  the  correction  factor  of 
the  measured  thermal  flux,  and  the  ratio  for  all  the  range  of  Rm  studied.  The 
same  relationship  is  valid  for  both  types  of  HFTs . 
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Figure  13:  Relationship  between  the  ratio  which  represents  the  collection  factor  of 

the  measured  thermal  flux,  and  the  ratio  RJRm,  for  the  physiological  range  of  Rm 
studied.  The  same  relationship  is  valid  for  both  types  of  HFTs. 


Tkble  3:  Tfemperature  and  heat  flux  data  for  the  last  15  minutes  of  a  3  hour  immersion  of  the 
forearm  in  water  at  different  temperatures.  Tw:  water  temperature  (°C);  Tn:  rectal 
temperature  of  the  subject  (°C);  Tc :  core  temperature  of  the  forearm  (°C);  Tsk'  uncovered 
skin  temperature  of  the  forearm  measured  with  a  fine  gauge  thermocouple  (°C);  Hwas: 
measured  heat  flux  from  the  skin  of  the  forearm  (W»m'2);  heat  flux  from  the  skin 
of  the  forearm  corrected  for  the  insulating  effect  of  the  HFT  (W»m‘2);  %  dev:  percentage 
of  deviation  between  Ameas  and 
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4.  DISCUSSION 

Because  of  their  convenience,  HFTs  are  being  used  increasingly  to  study  human  thermal 
physiology  under  different  environments!  conditions.  However,  the  vast  majority  of  studies  have 
used  the  factory  calibration  data  provided  with  the  HFTs,  and  little  attention  has  been  given  to  the 
effect  of  the  thermal  resistance  of  the  device  on  the  true  thermal  flux  from  the  skin.  The  results  of 
this  study  indicate  that  the  selection  and  the  use  of  HFTs  must  be  done  with  care. 


4.1  Calibration  of  HFTs 

The  recalibrations  of  HFTfc  from  two  companies  have  shown  different  deviation  values 
between  the  factory  supplied  and  the  experimentally  determined  calibration  constants.  The  results 
indicate  that  the  mean  difference  (in  absolute  value)  between  the  experimental  calibration  constants 
and  the  factory  constants  is  5  times  larger  in  the  case  of  Thermonetics  HFTs  (20.2%)  than  for 
Concept  HFlfc  (4.0%).  For  each  Thermonetics  HFT  calibrated,  the  factory  constant  was 
overestimated,  as  opposed  to  the  Concept  HFT  for  which  the  factory  constants  were  about  equally 
over  and  underestimated. 

A 

These  results  are  consistent  with  the  findings  of  Nuckols  and  Piantadosi  (7)  and  of  Snellen 
and  Frim  (11)  who  found  that  the  factory  calibration  constants  for  the  Thermonetics  HFTs  were  on 
average  9.0%  (3  to  20%)  and  26.1%  (8  to  40%)  overestimated,  respectively,  when  compared  to 
the  experimentally  determined  calibration  constants.  Sowood  (12)  found  that  the  factory 
calibration  constants  for  the  Concept  HFTs  differed  from  his  calibration  by  only  6  2%.  The 
advantage  of  the  present  study  over  the  others  is  that  the  calibrations  of  both  types  of  HFTs  were 
performed  by  the  same  reliable  and  reproducible  method  of  calibration  over  a  larger  range  of 
thermal  flux  (from  0  to  500  W*m'2),  making  the  comparison  between  the  two  types  of  HFTs  more 
accurate. 

The  differences  between  the  manufacturer's  and  the  experimental  values  may  be  due  to  an 
underestimation  of  the  heat  flux  by  the  manufacturer’s  method  because  the  transducers  were 
■  calibrated  prior  to  waterproofing  (9,12),  or  because  the  method  did  not  allow  for  the  insulative 
properties  of  the  disc  material  (12).  Whatever  the  cause  of  that  difference,  it  is  suggested  that  in 
the  case  of  HFTs  from  Concept  Engineering,  the  devices  can  be  used  without  recalibration,  since 
the  difference  between  the  manufacturer  and  the  experimental  calibration  data  is  small  (-5%). 
However,  potentially  large  errors  in  calculated  heat  flow  may  be  present  when  HFTs  from 
Thermonetics  are  not  recalibrated,  particularly  when  transducers  with  high-percentage  calibration 
errors  are  placed  at  heavily  weighted  skin  sites.  For  this  reason,  it  is  recommended  to  recalibrate 
individually  the  Thermonetics  transducers  before  using  them.  If  accurate  values  of  heat  flux  are 
required,  it  is  advisable  to  carefully  recalibrate  any  type  of  HFTs. 
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4 . 2  Correction  for  the  thermal  resistance  of  the  transducers 

The  study  involving  the  variable-R  model  indicates  clearly  that  a  strong  relationship  exists 
between  the  underlying  “tissue”  thermal  resistance  and  the  relative  error  in  thermal  flux  due  to  the 
thermal  resistance  of  the  HFT.  For  high  values  of  “tissue”  insulation  (Rt  >  0.5°C*m2»W'1;  Rf/Rm 
<  0.02;  see  Fig.  12),  the  measured  value  of  heat  flux  is  the  true  thermal  flux  from  the  tissue  (~0% 
error),  but  for  very  low  values  of  “tissue"  insulation  ( RT  <  2xl0‘4oC*m2»W‘1;  RjlRm  >  50;  see 
Fig.  12),  an  error  as  high  as  500%  can  be  present  in  the  measurement  of  the  thermal  flux  by  an 
HFT. 


lb  understand  this  effect,  take  the  example  of  two  systems  involving,  respectively,  a 
"nearly  perfect  conductor”  (i.e.  Rm<<RT)  and  a  "nearly  perfect  insulator"  (i.e.  Rm  >  >  Rf,  see 
Fig.  14).  If  an  HFT  is  fixed  on  the  conductor,  then  the  HFT  becomes  the  major  resistance  to  the 
heat  transfer  through  the  system  at  that  location.  Consequently,  heat  flow  measured  at  the  HFT  site 
(tineas) b®  much  less  than  through  the  conductor  only,  and  a  correction  is  necessary. 

On  the  other  hand,  when  a  "nearly  perfect  insulator"  is  used,  the  insulation  Rr  added  by  the 
HFT  to  the  system  (Rm  +  RT)  is  negligible  as  compared  to  the  insulation  value  of  the  insulator 
only.  Therefore,  the  heat  transfer  through  the  system  is  virtually  the  same  as  through  the  insulator 
only.  In  that  case,  the  HFT  will  closely  measure  the  true  value  of  heat  transfer  through  the 
insulator.  It  is  therefore  evident  that  the  difference  between  ti^^  and  ti^  is  highly  dependent 
upon  the  insulation  of  the  material  (or  tissue)  under  the  HFT.  The  deviation  from  the  true  value  of 
thermal  flux  increases  with  decreasing  insulation  of  the  underlying  material. 

The  above  interaction  becomes  particularly  important  when  HFTs  are  used  on  highly 
conductive  materials  such  as  metallic  skinned  mannequins  (copper  or  aluminum).  In  that  case, 
because  the  underlying  material  insulation  value  is  low  (RjfRm  >  1.5),  the  values  of  heat  flux  read 
by  the  HFT;  may  underestimate  the  true  value  of  heat  flux  by  more  than  100%,  depending  on  the 
insulation  of  the  mannequin  skin  (see  Fig.  12).  From  a  physiological  standpoint,  the 
underestimation  of  heat  flux  through  the  skin  measured  by  an  HFT  can  be  important  when  the 
device  is  used  during  vasodilatation  in  warm  subjects,  as  suggested  by  Gin  et  al.  (6)  and  Wissler 
and  Ketch  (14).  For  that  situation,  as  shown  in  Fig.  13,  the  underestimation  of  the  measured  heat 
flux  can  range  between  29  and  35%  from  the  true  value  of  heat  flux  (ratio  R  j!Rt  ranging  between 
0.45  and  0.6;  ref.  3),  depending  on  the  HFT  model  used  and  the  degree  of  vasodilatation  in  the 
tissue.  However,  when  HFTs  are  used  on  vasoconstricted  skin  in  cool  subjects  (ratio  RT/R( 
ranging  between  0.025  and  0.15;  ref.  3),  the  underestimation  of  the  measured  heat  flux  is  much 
less,  ranging  between  3  and  13%.  It  may  therefore  not  be  as  important  to  correct  the  heat  flux 
values  for  the  thermal  resistance  of  the  HFT  when  the  device  is  used  on  vasoconstricted  skin. 
However,  it  becomes  important  to  use  Wissler  and  Ketch’s  (14)  equation  to  correct  the  measured 


mees  corr 


Figure  14:  The  two  thermal  systems  used  to  explain  the  effect  of  the  underlying  tissue 
on  the  correction  factor  of  Hmas  (see  text  for  details). 
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heat  flux  for  the  thermal  resistance  of  the  HFT,  when  the  device  is  used  on  vasodilated  skin.  When 
the  correction  equations  are  used,  the  parameter  RT  in  the  equations  (thermal  insulation  of  HFT) 
has  the  value  of  0.01040C*m2*W'i  for  the  Thermonetics  HFTs,  and  0.0064°Om2»W*1  for  the 
Concept  Engineering  HFTs . 

Fig.  15  presents  a  practical  example  of  the  magnitude  of  the  error  due  to  the  thermal 
resistance  of  an  HFT  when  a  recalibrated  HFT  is  used  to  measure  the  thermal  flux  from  a  human 
forearm  immersed  in  water  at  different  water  temperatures.  The  heat  loss  from  the  forearm 
increases  with  the  reciprocal  of  water  temperature.  At  water  temperatures  below  the  lower  critical 
temperature  in  water  (~30°C;  see  ref,  10),  the  tissues  are  vasoconstricted  after  immersion  in  water 
(2).  The  %  error  due  to  the  thermal  resistance  of  the  HFT  is  only  9.8%  when  the  tissues  are 
vasoconstricted,  but  increases  to  25.3%  at  a  water  temperature  of  36°C  where  the  tissues  are 
partially  vasodilated.  These  data  are  in  agreement  with  the  results  found  with  the  variable-R 
model.  They  also  support  the  conclusion  that  the  value  of  thermal  flux  measured  by  an  HFT 
should  be  corrected  for  the  error  due  to  the  thermal  resistance  of  the  device  when  the  underlying 
tissues  are  in  a  vasodilated  state. 


4.3  Selection  and  use  of  HFTs 

Each  of  the  two  different  models  of  waterproof  HFTs  studied  presents  its  own  advantages. 
The  HFTs  from  Thermonetics  are  more  sensitive  than  those  from  Concept  Engineering.  This 
difference  might  be  noticeable  when  low  values  of  heat  flow  are  involved,  and  particularly  if  the 
voltage  outputs  from  the  HFTs  are  measured  by  a  voltmeter  with  low  sensitivity.  On  the  other 
hand,  the  HFTs  from  Concept  Engineering  don't  need  recalibration,  the  correction  for  the  thermal 
resistance  of  the  HFTs  is  less  important  (since  the  thermal  conductivity  is  almost  twice  that  of  the 
Thermonetics  HFTfe),  the  temperature  reading  from  the  embedded  thermistor  is  less  affected  by  the 
edge  effect  (see  ref.  14  for  details),  die  disc  is  less  susceptible  to  mechanical  breakdown.  But 
whatever  the  selection,  when  properly  used  and  calibrated,  both  HFTs  are  equivalent. 

lb  minimize  the  stress  applied  on  an  HFT  and  to  maximize  the  thermal  transfer  between  the 
skin  and  the  HFT,  the  device  should  be  fixed  properly  on  the  skin.  It  is  suggested  to  fix  the  HFT 
on  the  skin  with  thin  surgical  tape  with  the  wire  parallel  to  the  body  segment  in  the  proximal 
direction.  It  is  recommended  to  clean,  dry  and  wipe  the  skin  with  a  protective  dressing  (Skin  Prep, 
United  Medical,  FL)  prior  to  attaching  the  HFT.  It  is  important  that  the  complete  surface  of  the 
HFT  disc  makes  good  contact  with  the  skin.  An  "S"  is  made  with  the  wire  and  then  taped  to  the 
skin  to  minimize  the  stress  on  the  wire  (see  Fig.  16).  During  removal,  the  tape  is  pulled  from  the  ' 
side  of  the  HFT  (to  reduce  the  stress  at  the  junction  between  the  wire  and  the  HFT  disc)  in  the 
same  direction  as  the  orientation  of  the  hairs  on  the  skin.  To  avoid  an  air  or  water  gap  between  the 
HFT  and  the  skin,  and  to  maximize  the  thermal  transfer,  it  is  advised  to  apply  a  thin  layer  of 
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thermally  conductive  grease  to  the  surface  of  the  HFT  which  will  be  in  contact  with  the  skin. 


UJ./VO  M0|}'IB8H 


Skin 

Film  of  thermally 
conductive  grease 


Figure  16: 


Schematic  representation  of  the  suggested  method  to  fix  a  HFT  on  the  skin  to  minimize 
the  stress  on  the  HFT  and  optimize  the  heat  flux  from  the  skm  to  the  HFT. 
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5.  CONCLUSION 

The  calibrations  of  two  different  models  of  waterproof  HFTs  have  shown  significant 
differences  between  the  recalibration  and  the  manufacturers'  constants.  Furthermore,  a  significant 
difference  exists  in  the  error  of  calibration  between  the  two  manufacturers,  the  error  of  calibration 
(in  absolute  value)  being  5  times  larger  in  the  case  of  Thermonetics  HFTs.  Therefore,  it  is 
necessary  to  recalibrate  the  Ihermonetics  HFTfc  but  not  those  from  Concept  Engineering.  This 
becomes  an  important  criteria  for  the  selection  of  HFTs  for  someone  who  does  not  have  access  to 
calibration  facilities.  However,  when  properly  calibrated,  both  models  of  HFTs  are  equivalent. 
The  method  presented  in  this  study  to  calibrate  the  HFTs  using  the  Rapid-k  instrument  is  easy  to 
use,  reliable  and  reproducible. 

The  HFTs  inevitably  underestimate  the  thermal  flux  being  measured  because  of  the  effect  of 
the  thermal  resistance  of  the  device.  By  using  a  variable-R  model,  it  was  shown  that  the  magnitude 
of  the  underestimation  is  proportional  to  the  reciprocal  of  the  underlying  tissue  insulation.  The 
conditions  that  minimize  the  deviation  from  the  true  value  of  thermal  flux  are  a  low  perfusion  rate 
in  underlying  tissue  (vasoconstriction)  and  a  large  thermal  resistance  between  the  HFT  and  the 
„  environment.  The  ideal  condition  for  using  HFTs  is  during  the  measurement  of  the  thermal  flux 
from  a  dressed  subject  during  exposure  to  cold  air.  For  that  situation,  it  is  less  important  to  apply  a 
correction  for  the  thermal  resistance  of  the  transducers.  However,  the  use  of  HFTs  on  nude 
subjects  in  highly  conductive  environments  (such  as  wates,  or  high  pressure  air/or  helium-oxygen 
atmospheres),  and  particularly  under  conditions  that  do  not  produce  vasoconstriction,  may  cause 
an  underestimation  of  the  thermal  flux  by  the  transducer  of  30  to  35%.  For  those  particular  cases, 
it  is  important  to  correct  the  measured  values  of  thermal  flux  for  the  thermal  resistance  of  the  device 
by  using  Wissler  and  Ketch's  (14)  equation  of  correction.  The  heat  flux  data  from  experimental 
immersions  of  the  forearm  in  water  at  different  temperatures  confirmed  the  results  from  the 
variable-R  model:  the  underestimation  of  the  measured  heat  flux  is  ~10%  during  vasoconstriction 
and  approximately  ~30%  during  partial  vasodilation  of  the  tissues. 

The  HFT  is  a  device  which  represents  a  compromise  between  sensitivity  and  accuracy  of 
measurement.  If  used  properly,  it  provides  useful  information  within  an  acceptable  range  of 
accuracy.  The  error  in  the  measurement  of  the  thermal  flux  should  not  exceed  ~5%  when  care  is 
taken  with  respect  to  calibration  and  the  effect  of  the  thermal  resistance  of  the  transducer. 
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Program  used  for  calibration  of  HFTs 
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The  program  was  written  in  HP  Basic  Language  for  a  HP-85  computer.  The  commands 
used  in  the  program  are  for  the  following  peripherals:  a  HP-8290 1M  flexible  disc  drive,  a  HP- 
3495  A  scanner,  a  HP-3455A  digital  voltmeter,  a  59309A  HP- IB  digital  clock,  a  HP-2934A  printer 
and  a  HP-9872A  plotter. 


Ifl  1  FROG  NAME  HFTCAL 
20  !  pen  HFT  CALIBRATION 
30  !  BY  M  DUCHARME 
AO  ! .  THIS  PROGRAM  PRINTS.  PLOTS 
TEMP  t  HEAT  FLQL  DATA  FROM  P*Pi: 
-K  AND  --T 

i  :at«  ikiu:es  it  temps  &  2 

Hr'S 

!  aRCGRfit1  SCANS  CONTINUALLY  A 
ND  SUNS  DATA 

?0  I  1  MIN  SUMS  ARE  AUERAGEO  flUE 
R  N  SCANS  AND  PRINTED 
*  80  '  CHO-TH!  CH1®T,-  j  CH2«0raoid- 
k 

°0  !  CH6 :  eiriTHFT  1  :CH7:  emf  HFT2  J  CHS 
: emfHFTJ 

103  !  CH9:ei»fhFT4iCH20:TooHFTl:C 
H22: TdownHFTl :CH21 : Th 
110  I  Ch*3 i TuoHFT2  SCH24: TdownHFT 
2 

120  I  CH26 :  TuoHFT3 ;CH27: TdownHFT 
J 

130  I  CH29s  TuoHFT4 iCH30: TdownHFT 
4 

14  I  XlsMINIMUM  TIME  UALUE 
150  !  X2: MAXIMUM  TIME  UALUE 
160  I  YlsMINIMUM  UOLT  UALUE 
170  !  Y2SMAXIMUM  UOLT  UALUE 
180  !  ZlsMINIMUM  TEMP  IWLUE 
190  I  Z2:  MAX I MUM  TEMP  UALUE 
200  I  N* NUMBER  OF  SCAN/HIN 
210  I  R»NUM8ER  OF  MIN 
220  I  TjMONTH  VALUE  ON  CLOCK 
230  I  OJDAY  UALUE  ON  CLOCK 
240  I  HI: HOUR  UALUE  ON  CLOCK 
250  I  MiMINUTE  UALUE  ON  CLOCK 
260  !  SliSECONDE  UALUE  ON  CLOCK 
270  I  JjTIME  IN  SEC 
280  I  SiMEAN  Th  FOR  1  MIN 
290  I  UiMEAN  Tc  FOR  1  MIN 
300  I  UtMEAN  0  FOR  1  MIN 


310  I  UsUOLT  FOR  TEMP 

320  !  TCOiTENP  FOP  CHANNEL  C 

330  !  UliUOLT  FOR  0 

340  !  X3:PREUI0US  TIME 

350  I  X4: ACTUAL  TIME 

360  I  F:RATIO  U/(S-U) 

370  I  Nt:NAME  OF  TEST 
380  !  Ql: SCANNER  CHANNEL  ARRAY 
390  !  02: SCANNER  CHANNEL  ARRAY 
400  !  L: COUNTER  FOR  TEMP  LOOP 
410  I  B:COUNTER  FOR  UOLT  LOOP 
420  I  A: MEAN  UOLT  FROM  HFT1 
430  I  E:MEAN  UOLT  FROM  HFT2 
440  I  G:MEAN  UOLT  FROM  HFT3 
450  I  I:f€AN  UOLT  FROM  HFT4 
460  CLEAR 

470  LINPUT  •FILENAME- .Nl 
480  PRINTER  IS  701 
490  PRINT  -FILEJ-SNI 
500  PRINT 
510  PRINT 

520  PRINT  USING  740  :  "TIME*. “ME 
AN  Th" 

530  PRINT  USING  740  ;  -MEAN  Tc" . 
"MEAN  Q" 

540  PRINT  USING  750  :  "O/Th-Tc" 
550  PRINT  USING  740  ;  "VHFTr, 

h»  1  f 

560  PRINT  USING  740  ;  "VHFT2". 
570  PRINT  USING  740  J  "",*HFT3". 
580  PRINT  USING  740  f  "VHFT4". 
590  PRINT 

600  PRINT  USING  740  :  *VW 

610  PRINT  USING  750  :  «" 

620  PRINT  USING  740  f  "Tup".»Tdo 
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650  PRINT  USING  740  t  "TupVTdo 
«n"  ** 

640  PRINT  USING  740  I  "TupVTdo 
un".W 

650  PRINT  USING  740  I  "TupVTdo 
660  PRINT 

670  PRINT  USING  740  S  "(«in)"."( 
deg  C)" 

680  PRINT  USING  740  ;  “(deg  C)», 
"(W/*2)" 

690  PRINT  USING  750  !  "(U/*2  deo 
0“ 

700  PRINT  USING  740  s  "(dog  O". 
"(dog  C)V<»V)" 

710  PRINT  USING  740  i  "(doo  C)V 
"(deg  C)","(mU)“ 

720  PRINT  USING  740  i  "(deg  C)*f 
"(dog  OVUM"  « 

750  PRINT  USING  740  i  “(doo  Cl", 
"(dog  C)","(i»U)" 

740  IMAGE  i,17(12A) 

750  IMAGE  f,2(17A) 

'  760  PRINT 
770  PRINT 
780  PLOTTER  IS  705 
790  LOCATE 

800  Xl-0  !  MINIMUM  TIME  UALUE 
010  X2-500  !  MAXIMUM  TIME  UALUE 
820  Yl-0  I  MINIMUM  0  UALUE 
850  Y2-450  I  MAXIMUM  0  VALUE 
840  Zl-10  I  MINIMUM  TEMP  UALUE 
850  Z2-IO0  I  MAXIMUM  TEMP  UALUE 
860  I  »«  X  AXIS  LABEL  (TIKE) 

870  PEN  1 

880  SCALE  XI, X2.Y1,Y2 
890  AXES  10. 10, XI. Yl, 6, 5 
900  LORG  6 
9)0  CSIZE  4 

*  920  FOR  X-Xl  TO  X2  STEP  20 
950  HOUE  XfYl-(Y2-Yl)*,0i 
940  LABEL  X 
950  NEXT  X 

960  I  Y  AXIS  LABEL  (0) 

970  LORG  8 

980  FOR  Y-Yi  TO  Y2  STEP  50 
990  HOUE  XI ,Y 
1000  LABEL  Y 
1010  NEXT  Y 
1020  LORG  4 

1050  HOUE  X1MX2-X1)/2.Y2»(Y2-Y1 

)".01 


1040  LABEL  Nt 

1050  I  y  axis  label  (tempi 

1060  PEN  2 

1070  SCALE  X1.X2.Z1.Z2 
1080  AXES  X2fl,X2,Zl,1.5 
1093  LORG  2 

1100  FOR  Y-Zl  TO  Z2  STEP  5 

1110  MOVE  X2.Y 

1120  LABEL  Y 

1150  NEXT  Y 

1140  PEN  1 

1150  DEG 

1160  LOIR  90 

1170  LORG  5 

1180  MOVE  Xl-(X2-X1)*.0?.Z1MZ2- 
Zl)/2 

1190  LABEL  "O  (H/»2)" 

1200  LD1R  0 
1210  LORG  4 

1220  MOVE  Xl+(X2-Xl)/2.Zl-(22-Zl 
)M 

1250  LABEL  "TIME  (i»in)“ 

1240  PEN  2 
1250  LOIR  90 
1260  LQR6  5 

1270  MOVE  X2*(X2-X1)».07.Z1*(Z2- 
Zi)/2 

1280  LABEL  "TEMPERATURE  (deg  0" 

1290  PEN  0 
1500  CLEAR 

1510  DISP  "PRESS  'CONTINUE'  TO  S 
TART" 

1520  PAUSE 
1550  I 
1540  I 
1550  R-0 
1560  N»0 

1570  OUTPUT  722  ;"F1R7TJM0H1" 

1580  CLEAR 

1590  OUTPUT  708  ("PR* 

1400  OUTPUT  708  J"T» 

1410  R-R*l 

1420  GOSUB  1560  I  SCANSTART 
1450  GOSUB  1660  !  SCAN 
1440  ENTER  708  \  T.O.H1.M.S1 
1450  J«SU6O"HO60O*H1 
1460  IF  J/60<R  THEN  1450 
1470  GOSUB  1970  I  CALSCAN 
1480  GOSU8  2060  !  PRINT 
1490  IF  R>1  AND  R  MX  5-0  THEN  G 
OSUB  2190  !  PLOT 


49 


1500  GOSUB  2540  UPDATE 

1510  GOTO  1410 

1520  STOP 

1550  i 

1540  ! 

1550  !  SCANSTART 
1560  S>0 
1570  U-0 
1580  U-0 
1590  N-0 
1600  A*0 
1610  E-0 
1620  G*C 
1630  I»0 
1640  RETURN 
1650  I  SCAN 
1660  N»N*1 
1670  !  TEHP  .SCAN 
1680  DIM  01(6) 

1690  DATA  21,20,22,23,24,26.27.2 
9,30 

1700  NAT  READ  01 

1710  FOR  L-0  TO  8 

1720  OUTPUT  709  USING  ‘'ODD'  !  Cl 

(L) 

1730  OUTPUT  722  {W 

1740  ENTER  722  :  U 

1750  T(LWU»1000)*2*(-.6?1)*25, 

897*1000*U-.01 

1760  NEXT  L 

1770  S*S*T(0) 

1780  IMJ+(T(l)*T(3)*T($)*T(7))/4 

1790  I  UOLT.SCAN 
1800  DIN  02(4) 

1810  DATA  2, 6, 7, 0,9 
1820  NAT  REAO  02 
1830  FOR  B-0  TC  4 
1840  OUTPUT  709  USING  "ODD*  j  02 
(8) 

1850  OUTPUT  722  )"UT5‘ 

1860  ENTER  722  ;  U1 
1870  u(B)-umooo 
1800  NEXT  B 
1890  U-Utg(0) 

1900  A»A*U(1) 

1910  E-E+UC2) 

1920  G*G*U(3) 

1930  I-I«V(4) 

1940  RESTORE 
1950  RETURN 
1960  !  CALSCAN 


WO  S-S/N 
1980  U«U/N 

1990  td«- .  0044*2? .  5  092*  ( U/N ) 

2000  A»A4) 

2010  E-E/N 
2020  G-G/N 
2030  IM/N 
2040  RETURN 
2050  !  PRINT 
2060  PRINTER  IS  701 
2070  F-U/IS-U) 

2080  PRINT  USING  2130  5  R.S.U,W, 
F 

2090  PRINT  USING  2140  i  T(1).T(2 
).A 

210C  PRINT  USING  2140  :  T(J).T(4 
).E 

2110  PRINT  USING  2140  i  T(5).T(6 
).G 

2120  PRINT  USING  2140  !  T(7).T(8 

).I 

2130  INAGE  *.3D.D.7X,3(3D.3D.5X> 
.3D.3D.10X 

2140  INAGE  ♦,3(3D.3D,5X) 

2150  PRINT 
2160  RETURN 
2170  PRINT 
2ieo  !  PLOT 
2190  X3-R-5 
2200  X4-R 

2210  SCALE  X1,X2,Y1,Y2 
2220  HOUE  X4,Y1 
2230  ! 

2240  PEN  4 
2250  HOUE  X3,U0 
2260  DRAW  X4.U 
2270  PENUP 

2280  IF  X4  NOO  30>5  OR  X4  HOD  30 

<5  THEN  2330 

2290  LOIR  0 

2300  LORG  5 

2310  HOUE  X4,U 

2320  LABEL  "O* 

2330  I 

2340  SCALE  X1,X2,21,Z2 
2350  PEN  2 
2360  HOUE  X3,S0 
2370  ORAU  X4,S 
2380  PENUP 

2390  IF  X4  NOO  30<15  OR  X4  NOO  3 
0>15  THEN  2420 
2400  NOUE  X4.S 

* 
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r 

m  2410  LABEL  “Th" 

r-  2420  I 

24 JO  PEN  J 
2440  tm  XJ.U6 
;  2450  DRAW  X4,U 

2460  PENUP 

%  2470  IF  X4  MOO  J0<25  OR  X4  MOO  3 

0>25  THEN  2500 
2480  MOUE  X4.U 
2490  LABEL  "Tc" 

2500  1 
2510  PEN  0 

•  2520  HOUE  X2,21 
25J0  RETURN 
2540  I  UPDATE 
2550  N-0 

2560  IF  R-l  OR  R  NOD  5-0  THEN  25 
80 

#  2570  RETURN 
2580  U0-U 
2590  UO-U 
2600  SO-S 
2610  RETURN 
2620  END 
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ANNEX  B 

Equivalence  of  propoced  equations  for  the  correction 
for  the  thermal  resistance  of  HFT 


Fig.  B-l  presents  the  system  involved  in  the  demonstration  of  the  theoretical 
equivalence  of  the  equations  proposed  by  different  authors  to  correct  for  the  thermal 
resistance  of  HFTs.  The  system  represents  the  immersion  of  a  human  body  segment  in 
water,  but  it  can  also  be  applied  in  air  with  consideration  for  the  boundary  layer.  In  this 
system,  l)corr  is  the  thermal  flux  (W  •  nr2)  through  the  skin  corrected  for  the  insulating 
effect  of  the  HFT,  which  corresponds  to  the  true  thermal  flux  through  the  skin  not  covered 
by  a  HFT.  ii^as  *s  the  measured  heat  flow  value  (W  •  nr2),  Tc  is  the  core  temperature  of 
the  segment  for  which  the  heat  flow  is  measured  (°C),  is  the  temperature  of  the 
uncovered  skin  of  the  segment  (°C),  7/  is  the  temperature  of  the  skin  under  the  lower 
surface  of  the  HFT  (°C),  Tu  is  the  temperature  of  the  upper  surface  of  the  HFT  (°C),  Tw  is 
the  temperature  of  the  water  (°C),  Rf  is  the  thermal  resistance  of  the  HFT  (°C  •  m2  •  W*1), 
R,  is  the  thermal  resistance  of  the  tissue  (°C  •  m2  •  W*1),  and  Rw  is  the  thermal  resistance  of 
the  water  (°C  •  m2  •  W-1). 


Weter 


HFT 

Skin 

Tissue 


Figure  B-l .  Schematic  representation  of  the  thermal  system  involved  in  the  demonstration 
of  the  theoritkal  equivalence  of  the  equations  proposed  by  different  authors  for  the 
correction  for  the  thermal  resistance  of  the  HFTs. 


From  Fourier's  law  of  heat  transfer  by  conduction: 


H 


Tc  ’  Tsk 


COIT  R, 


H 


^sk  *  Tw 


coir  R 


w 


H 


cotr 


Tc-Tw 
Rt  -  Rw 


In  addition,  Anucs  can  be  described  by: 


H 


Tc  *  Tj 


meas  r 


H 


Ti  -  T„ 


meas  R. 


H 


T  -  T 
u  w 


meas  R 


w 


H 


T  -  T 
1  c  w 


meas  Rt  +  RT  +  Rw 


Knowing  that  from  equation  B-5: 


Tu  “  Tj  -  RT  •  H 


meas 


(B-l) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 


and  from  equation  B-6: 


T  -T 

„  u  w 

R  „  — ; - 

W  H 

meas 


(B-9) 


we  can  now  express  Eq.  B-2  as  follows: 


H 


^Tsk  ’  V  *  Hmeas 


(Tj-Tw)-Rt.H 


meas 


(B-10) 


Eq.  B-10  was  proposed  by  Strong  et  al.  (1985). 
From  Eq.  B-4: 


V 

H 


meas 


(B-ll) 


and  from  Eq.  B-5: 


Tj  =  Tu  +  R7  •  Hmeas 


(B-12) 


Therefore,  Eq.  B-l  can  now  be  expressed  as  follows: 


H 


<V  *  ^meas 


COIT 


(B-l  3) 


(Tc-T,) 
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H 


Tc '  Tsk 


00,1  (Tc-  Tu)/Hmcas  "  Rt 


(B-14) 


Eqs.  B-13  and  B-14  are  simply  rearranged  forms  of  the  equation  of  Gin  et  al.  (1980). 
From  Eqs.  B-3  and  B-7: 


&corr  *  (Rs +  Rw)  °  Hmeas  *  (Rt +  RT  +  Rw> 


(B-15) 


^corr  *  (Rt +  Rw)  *  fyneas  *  (Rt +  Rw)  +  fyneas  *  RT  (B- 16) 


By  rearrangement,  Eq.  B-16  becomes: 


^rneas  '(Rt+Rw) 


H«»r,(Rt+Rw>*Hmeas  -RT 


(B-17) 


H 

con- 


**mcas  *  ^cotr  *  (Rt+Rw) 
Hcon-  *  (Rt  +  Rw)  'Hmcas  •  Rt 


(B-18) 


By  using  Eq.  B-3,  Eq.  B-18  can  be  rearranged  and  expressed  as  follows: 


H 


H 


meas 


corr 


1  ’^mcas  *  rt/(Tc"  Tw^ 


(B-19) 


Eq.  B-19  is  the  correction  equation  derived  by  Wissler  and  Ketch  (1982),  and  later  by 
Sowood  (1986). 
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This  mathematical  demonstration  proves  the  theoritical  equivalence  of  all  the  equations 
proposed  by  different  authors  to  correct  for  the  thermal  resistance  of  riFTs. 
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